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We investigate dynamical properties of bright solitons with a finite background in the 
F = 1 spinor Bose-Einstein condensate (BEC), based on an integrable spinor model which 
is equivalent to the matrix nonlinear Schrodinger equation with a self-focusing nonlineality. 
We apply the inverse scattering method formulated for nonvanishing boundary conditions. 
The resulting soliton solutions can be regarded as a generalization of those under vanishing 
boundary conditions. One-soliton solutions are derived in an explicit manner. According to 
the behaviors at the infinity, they are classified into two kinds, domain-wall (DW) type and 
phase-shift (PS) type. The DW-type implies the ferromagnetic state with nonzero total spin 
and the PS-type implies the polar state, where the total spin amounts to zero. We also 
discuss two-soliton collisions. In particular, the spin-mixing phenomenon is confirmed in a 
collision involving the DW-type. The results are consistent with those of the previous studies 
for bright solitons under vanishing boundary conditions and dark solitons. As a result, we 
establish the robustness and the usefulness of the multiple matter- wave solitons in the spinor 
BECs. 

KEYWORDS: spinor Bose-Einstein condensate, internal degrees of freedom, matrix nonlinear 
Schrodinger equation, integrable model, bright soliton, nonvanishing boundary 
condition, multi-component Gross-Pitaevskii equation, spin exchange coupling, 
spin switching 



*E-mail address: kurosaki@monet.phys.s.u-tokyo. ac.jp 



1/18 



J. Phys. Soc. Jpn. Full Paper 

1. Introduction 

In 2002, matter-wave bright solitons in quasi-one-dimensional (ID) Bose-Einstein conden- 
sates (BECs) were observed experimentally.^'^) Bright solitons propagate in most eases with 
much larger amplitudes than dark solitons,^' and are expected to have the potential for 
various applications such as coherent transport and atom interfcronietry. Soliton propagation 
in BEC can be described by the Gross-Pitaevskii (GP) equation. The CP equation, called the 
nonlinear Schrodinger (NLS) equation in nonlinear science, is integrable and has soliton solu- 
tions in a one-dimensional and uniform system. Recent experimental and theoretical advances 
about matter-wave bright solitons are reviewed, for instance, in ref. 5. 

The experimental creation of matter-wave solitons has been so far achieved only for single- 
component BEC. It is, nevertheless, very interesting to consider soliton propagation in BEC 
with internal degrees of freedom, so-called, spinor BEC. When BEC of ultracold alkali atoms 
is trapped exclusively by optical means, the hyperfine spin of atoms remains liberated. The 
spinor BEC was realized in such a way.^"^) Internal degrees of freedom endow solitons with 
a multiplicity. The multiple solitons will show a rich variety of dynamics. Here, we focus on 
the boson system in the F = 1 hyperfine spin state, exemplified by ^'^Na, ^^K and ^'^Rb. The 
multi-component GP equation for F = 1 spinor BEC turns to an integrable model at special 
points, which is mathematically equivalent to the matrix NLS equation. An integrable model 
with a self-focusing nonlineality enables one to perform exact analysis via the inverse scatter- 
ing method (ISM) for the matrix NLS equation.^) In particular, bright soliton solutions under 
vanishing boundary conditions (VBC) are obtained, whose properties are investigated in refs. 
10 and 11. Recently, the ISM for the matrix NLS equation under nonvanishing boundary con- 
ditions (NVBC) is formulated. ""^^^ Dark solitons in the F = 1 spinor BEC can be investigated 
by applying the ISM under NVBC to an integrable model with a self-defocusing nonlineal- 
ity. ""^^^ Although the ISM under NVBC is dedicated mainly to the self-defocusing case, we note 
that this technique is also applicable to an integrable model with a self-focusing nonlineality, 
which makes us available to bright soliton solutions with a finite background. 

In this paper, the detail of matter-wave bright solitons in the quasi-lD F = 1 spinor BEC 
is further investigated, based on an integrable model. We consider matter-wave spinor bright 
solitons traveling on a finite background of the condensate. We write down explicitly new 
soliton solutions, and verify that the obtained soliton solutions have the similar properties 
compared to those without a background. In the usual experimental setups, the condensates 
are confined in a finite-size regime, and the matter-wave bright solitons will accompany a 
finite background. The study given in this paper is meaningful in such realistic circumstances. 

The paper is organized as follows. In § 2, the GP equation for quasi-lD F = 1 spinor BEC 
is introduced. In particular, the integrable model is presented. There, the interactions between 
two atoms are supposed to be inter-atomic attractive and ferromagnetic, which lead to bright 
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solitons. In § 3, the inverse scattering method under nonvanishing boundary conditions is 
apphed to the integrable model. This apphcation leads to bright soliton solutions with a finite 
background. Several conserved quantities of the model are also provided. One-soliton solutions 
are investigated in § 4. The spin states of one-solitons are classified, assuming that discrete 
eigenvalues are purely imaginary. Two-soliton solutions are discussed in § 5. The last section, 
§ 6, is devoted to the concluding remarks. 

2. GP Equation for F — 1 Spinor Bose-Einstein Condensates 

For BEC of ultracold alkali atoms, the mean-field theory works well, because almost all 
atoms go into condensation and the condensate is dilute. In this paper, we deal with the 
quasi-one-dimensional system. Atoms in the F = 1 hyperfine spin state have three magnetic 
substates labeled by the magnetic quantum number mp = 1,0,-1. The system is charac- 
terized by a vectorial field operator with the components corresponding to each substate, 
^ = ($1, $0) ^-i)"^) satisfying equal-time commutation relations: 

where the subscripts a, (3 take on 1, 0, ~1. In the framework of the mean-field theory for BEC, 
the quantum field is replaced with the order parameter: 

^{x,t) is often called the spinor condensate wavefunction, which is normalized to the total 
number of atoms iVr: 

/dx*t,.,,.„.,, = A,,. (3) 

The spinor condensate wavefunction obeys a set of coupled evolution equations, namely, the 
multi-component GP equation: 

in^t^i = -^dl^i + (co + C2) (|$i|' + l^ol") $1 

+ {co - C2)|$-lP$l + C2$li$0' 

ihdt^o = -^9l^o + (co + C2) + l^-iP) ^0 

-Fco|$oP*o + 2C2$0^1$_1, 
ihdt^.l = -^dl^.i + (co + C2) + |$oP) $-1 

+ (co-C2)|$lP$-l + C2$t$0> (4) 

where cq = {go + 2^2 )/3 and C2 = (52 — 5o)/3 denote effective ID coupling constants for the 
mean-field and the spin-exchange interaction, respectively. Here, the effective ID coupling 
constants gf are given by-*^^) 

_ _ 4^^a/ 1 , . 

mall-Caf/a^' 
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where a/ are the ,s-wave scattering lengths in the total hyperfine spin / channel, a± is the 
size of the transverse ground state, m is the atomic mass, and C = —(^(1/2) ~ 1.46. Note that 
one may change the values of cq and C2 by tuning a±. 

Equation (4) is derived as follows. The interaction between two atoms in the F = 1 
hyperfine spin state has a form,^^'^^) 

V{X1 - X2) = S{XI - X2) (co + C2F1 • F2) , (6) 

where Fj is the spin operator. The Gross-Pitaevskii energy functional is thus given by 

Eopm = J dx(£^d,Kd.^a + ^KK'^a'^a + ^KK'^afi ' ^a' 0' ^ (3' ^ (3^ , (7) 

where repeated subscripts {a,P,a',P' = 1,0,-1) should be summed up and f = (f^,f^,f^)-^ 
with P{i = x,y,z) being 3x3 spin-1 matrices. Then, the variational principle: ihdt^o:{x,t) = 
dEGp[^]/S^l{x,t), for a = 1,0,-1, yields eq. (4). 

An important fact is that eq. (4) possesses a completely integrable point when cq = C2 = 
— c < 0, equivalently, 2go = —§2 > 0.^^'^-*^^ This condition is realized when 

a± = 3C- — - — , (8) 
2ao + 02 

assuming that 0002(^2-00) > holds. The situation corresponds to attractive (cq < 0) and fer- 
romagnetic (c2 < 0) interaction. When we change the wavefunction by $ = (^1, \/2^0) 0-i)"^ 



and measure time and length in units of t = ha±/c and x = fi^ a^j^mc, respectively, we can 
rewrite eq. (4) with cq = C2 = — c < into the dimensionless form, 

15*01 = -dl<^i - 2(|</'iP + 2|(/.op)</.i - 2d'^g, 

\dt<\>^ = -dlcPo - 2(|</.-ip + l^ol' + \(pi\^)<Po - 2(t>lM-u 

idt4>-i = -dl4>-i - 2(|</)_i|2 + 2|(/)o|2)(/,_i - 24>14>1. (9) 

Then, eq. (9) is found to be equivalent to a 2 x 2 matrix version of the NLS equation with a 
self-focusing nonlinearity: 

idtQ + dlQ + 2QQtQ = o, (10) 

with an identification, 

V <Po 9-1 J 

The matrix NLS equation (10) is integrable in the sense that the initial value problem can 
be solved via the inverse scattering method.^' The integrability of the reduced equations 
(9) is thus proved automatically. Thus, we have derived the integrable spinor model. Another 
integrable point of eq. (4) is cq = C2 = c > 0, i.e., the matrix NLS equation with a self- 
defocusing nonlineality.^^^ Special solutions for generic coupling constants cq, C2 are given in 
ref. 17. 
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Wc consider bright sohton solutions of tlic intcgrablc model (9) under NVBC, whereas 
those under VBC are studied in refs. 10 and 11. We summarize briefly the results of the inverse 
scattering method for eq. (9) with NVBC.^^^ 

We define the nonvanishing boundary conditions as 

Q{x,t) —>■ Q±, X —>■ ±oo, 

QiQ± = Q±Qi = xp, (12) 

where Aq is a positive real constant and / denotes a 2 x 2 unit matrix. Note that vanishing 
boundary conditions are recovered as Aq 0. The analysis of the ISM under NVBC yields 
the standard form of the multiple soliton solutions of the 2x2 matrix NLS equation with a 
self- focusing nonlineality (10) as 

/ Hie^i \ 

I + 2i{I---I)S-^ 



Q{x,t)=Xoe"^{ 



N 



(13) 



Here, a 2 x 2 complex matrix Ilj is called the polarization matrix. S is a 2N x 2N matrix 
defined by 

Ao 



Cj + 



SijI + 7 



\2 



i(Ci + 0) VCi + A 



1 



+ 



1 



n,:e>^% 



+ A,-. 

l<i,j<N, 



(14) 



where Xj is a complex discrete eigenvalue for the bound state and = {X'j + Aq)^/^ with 
Im Q > for j = 1, . . . , A^. It is required for the ISM under NVBC that a two-sheet Riemann 
surface is introduced appropriately, due to a double- valued function Q. The phase of the 
carrier wave, (f){x,t), is given by 

(f){x, t) = kx- {k^ - 2Xl)t + 6, (15) 

and the coordinate function is given by 

Xj = Xji^, t) = 2iO(x - 2(A, + k)t). (16) 

The above solution is the M{= iV/2)-soliton solution. The ISM under NVBC for the self- 
focusing case results in pairs of discrete eigenvalues corresponding to each Riemann sheet. 
The constraint should be imposed on Xj and Q [j = 1, ■ ■ ■ ,N) such that X21-1 = X21 and 
C,2l-i = ~C2i / = !,•■■ , N/2. At the same time, II^' must satisfy that ^21-1 = 112^. 

For our reduction to the integrable model for F = 1 spinor BEC, we must make the 
potential Q symmetric, noting eq. (11). The symmetry of Q is naturally reflected in Hj. 
When we take every Yij to be symmetric in eq. (13), soliton solutions of the integrable model 
(9) under NVBC are obtained. 
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The form (13) is the standard form of sohton solutions in the sense that the boundary 
value at a; — > oo -^^^ e^^ -^^ is supposed to be fixed as 

X — oo. (17) 

The spinor model, however, allows the SU{2) transformation of the solutions, if they are kept 
symmetric. To be concrete, let W be a 2 x 2 unitary matrix. When Q is a solution of eq. (10) 
with eq. (11), then 

Q' = UQU^, (18) 

is also a solution. Assuming that Q is the standard form (13), the limit x —>■ oo oi Q' becomes 

Q' e-^-^^^'*) ^ Ao WW^ = \oQ'+, X ^ oo. (19) 

Q'j^ = UU^ is the so-called Cholesky decomposition. The arbitrary boundary conditions Q'j^ 
other than eq. (17) are thus realized via the SU{2) transformation. On the other hand, the 
behavior in the limit x ^ —oo varies depending on whether detll = or not, which will be 
discussed later for the one-soliton case. 

There is another important concept about an integrable model. Due to the integrability, 
the model has the infinite conservation laws, which restrict the dynamics of the system in an 
essential way. Several conserved quantities, related to the physical quantities of the system, 
are listed below: 

Total number : iVr = J dx n{x,t), 

n(x,t) =tr(Qtg) _tr(Q^Q±). (20) 
Total spin : Fx = J dx f{x,t), 

/(x,t) =tr(Qto-Q). (21) 
Total momentum : Px = J '^^ Pi^^'t)-, 

p{x, t) = -in[tr (QtQ^) _ tr {QiQ±,^)]. (22) 
Total energy : Et = J dx e{x,t), 

e{x,t) = c[tr {QIq^ - Q^QQ^'Q) 

-tr {Qi,^Q±,^ - QiQ±QiQ±)]. (23) 

Here, cr = [a^, (Jy, (Jz)'^ are the Pauli matrices. To avoid the divergence of integrals, wc should 
subtract the contribution of the background from the physical quantities, except the total spin, 
to which the background does not contribute explicitly. These subtractions arc emphasized 
by the bars on the conserved densities and quantities. The local spin density / = (f^, fy, fzY' 
is covariant under the SU{2) transformation (18), whereas the other densities such as n(x,t), 
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p{x, t) and e{x, t) arc invariant. The total macroscopic spin is directed to face to an arbitrary 
direction by the global spin rotation. The SU (2) symmetry of the system causes the energy 
degenerated states of solitons for this spin rotation. 

4. One-Soliton Solutions 

In this section, one-soliton solutions of the integrable spinor model (9) are investigated in 
detail. We can derive the explicit form of one-soliton solutions by setting N = 2 (M = 1) in 
the formula (13). The calculation is complicated but straightforward. The result is as follows: 

T 



(24) 



g = Aoe^^(-*)(/ + 2.^^^^ 
where det S is given by 

detS = k^kI — e-^ViKiK^ tr n — e^^f2KiK2 tr 11^ 
+ eXi+X2^^^2 {tu + i/iz/2(|trn|^ - 1)} 
+ e^^i i^fnl det U + e^^Vl^f det W 
- e^^i+^2j^iK2^^ tr det U - 6^^+"^^^ iy2Kiw tr U det 
^g2xi +2x2^2 i^g^np, (25) 
and T is a 2 X 2 matrix such that 

_exi+x2^^^2|^itrn-nt + ^2trnt -n + z^dtrnp- 1)/| 

- c^xi /iiK^ det n ■ / - c^^^ ^121^1 det H'f • / 
+ e2»+>^2^2 detnj^i^n^ +tutrn^ •/} 
+ e^i+2^2^1 det nt {^1 n + TO tr n • /} 

_g2xi+2x2^(^^ + ^2)|detnp.7. (26) 

We explain physical meanings of notations. The phase of the carrier wave is given by 

(j){x, t) = kx- {k^ - 2\l)t + 5. (27) 

Let \j and (j = (A^ + '^o)^'^^ 3 — be complex constants satisfying Ai = A^ and 
Ci = —CI- Without loss of generality, we assume that (Ai,Ci) ((A2,C2)) belongs to the upper 
(lower) Riemann sheet, which is characterized such that Im ^ Im A > (Im ^ Im A < 0). 
Xj{x,t) is expressed in terms of them as 

Xi = Xi{x, t) = 2iCi(x - 2(Ai + fc)t), (28) 

X2 = X2{x, t) = 2iC2(x - 2(A2 + k)t). (29) 

Note that xi = X2 = X holds. Rex thus denotes the coordinate of the envelope soliton, 
whereas Im x implies the self-modulation phase. The polarization matrix 11 is a 2 x 2 symmetric 
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matrix. Here, the normalization in a sense of the square norm is imposed on 11 as a matter of 



convenience: 



n^(^"V 2|ap + |^p + |7p = l. (30) 
V " 7 / 

The other parameters are expedient functions of Aq, Ai and A2: 

Kj = ———, ll = l\o , Vj = —^, W = ViV2-H, <ij=Uj-H, (31) 

Ci + Aj Ci + C2 

for j = 1,2. We Ust the meaning of each parameter as follows: 

k : wave number of soliton's carrier wave. 

Ao : amplitude of soliton's carrier wave. 

(f){x, t) : phase of soliton's carrier wave. 

Re x{x,t) : coordinate of soliton's envelope. 

Im x{x,t) : self-modulation phase of soliton. 

n : symmetric polarization matrix of soliton. 

Equations (24)- (26) are new soliton solutions which had never been written down explicitly 
in the literatures. If we take the vanishing limit Aq — > 0, and (2 converge at Ai and — A2, 
respectively. Then, 

O (i^ eA . 2k. ne-(x.+^/^) + (a^nta^)eXH+p/2detn 

VlAil ; e-(2x«+P) + i+e2x«+P|detn|2 ^ ' ^^^^ 



with notations 



(33) 



2kR' 

XR = XR{x,t) = kR{x - 2kit) - e, (34) 

XI = Xi{x, t) = kjx + {k\ - kf)t, (35) 

each of which holds the following correspondence respectively: 

kn = -2ImAi, (36) 

ki = 2ReAi + k, (37) 

e = -ln(4|Ai|). (38) 

Equations (32)- (35) are the same forms as those in ref. 11, except a phase factor. This con- 
sequence is natural but non-trivial, because the formula of solitons under VBC^^ is quite 
different from that under NVBC (13), in particular, in the form of the matrix S. Actually, the 
initial displacement e can be arbitrarily changed, regardless of eq. (38), by the parallel shift 
of the position x. We have shown that the soliton solutions (24)-(26) can be regarded as a 
general form of bright soliton solutions, including the case of VBC. 
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4.1 Classification by the boundary conditions 

Wc shall show that there arc two kinds of one-soliton solutions depending on the boundary 
conditions, dct 11 = or det 11 7^ 0. The similar classification about the boundary conditions 
also exists for dark solitons.-*^'^) The examples of snapshots of one-soliton density profiles are 
shown in Fig. 1. The upper row is for detH = 0, and the lower row is for dctll 7^ 0. The 
shape of envelope solitons looks a locally-oscilating wave rather than, literally, a solitary wave, 
because of the self-modulation due to the complex velocity. 

For detn = 0, the boundary conditions of the standard form (24)-(26) are 

Q e""^ — Aq/, X ^ 00, 

?i tr n • nt + <;2 tr nt ■ n + /x(|tr np - 1)7" 



ge-^<^^ Ao 



2i 



-00. (39) 



W + UiU2{\tTU\'^ - 1) 

The left and right boundary values differ in not only the global phase but also the population 
of each component, in general. That is, those are the SU{2) rotated boundary conditions. 
In the upper row of Fig. 1, we see that the envelope soliton of each component forms the 
domain-wall (DW) shape, although it does not manifest in the total number density. 
On the other hand, for det 11 7^ 0, the boundary conditions are 

Q e""^ —>■ XqI, X —>■ 00, 

Q e-'"^ ^Xq(i- 2i ^^^^ I, -00. (40) 

In contrast to the case that det 11 = 0, both boundary values are diagonal matrices, and only 
the phase-shift (PS) occurs. That is, those are the U{1) rotated boundary conditions. 

For the above reasons, we call one-soliton solutions the DW-type for det 11 = 0, and the 
PS-type for det 11 7^ 0. Remark the following; the spin density profile of DW-type suggests 
that the total spin is nonzero, whereas that of PS-type is dipole-shape, implying that the total 
spin amounts to zero. See the right panel of Fig. 1. This observation will be solidified in § 4.3. 
4-2 Case of purely imaginary discrete eigenvalues 

The ISM performed on Riemann sheets involves a double-valued function of the spectral 
parameter, and it usually renders a very complicated representation of 7V-soliton solutions 
even for N = 1, ss is seen from eqs. (24)-(26). To simplify an explicit representation, it is 
convenient to assume that A^ and C,j are purely imaginary. The similar approach is employed 
for the analysis of A/"-soliton solutions of the derivative NLS equation under NVBC.^°) 

If we take a pair of discrete eigenvalues as 

(Ai, Ci) = (iAoA, iAoC), (A2, C2) = (-iAoA, iAoC), (41) 
where A and C are positive real numbers such that 

A > 1, C = - 1, (42) 
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Fig. 1. Snapshots of onc-soliton density profiles. The upper row is plotted for det 11 = at the moment 
t = 0, with = 0, Ao = 1, Ai = 1+i, = 1.27 + 0.791 ix{x,t) = -(1.57- 2.54i)a;- (8.23- 1.94i)t) 

'4/5 2/5' 



and n 
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The lower row is plotted for det 11 ^ at the moment t = 0, with the 



/ I/V2 2/5 \ 

same parameters except for 11 = _ . The left panel (a) depicts the local density 

for each component, \<pi\'^ (solid line), \(f)o\'^ (chain line) and |<;^-i|^ (dotted line). The center panel 
(b) depicts the local number density n, where the contribution of the background is included. The 
right panel (c) depicts the local spin densities, fx (solid line) and fz (dotted line), fy vanishes 
identically due to a choice of a real matrix 11. 



we obtain a relatively simple form of one-soliton solutions as 
g^Aoe^^(^.*)(/ + 2i^), 



(43) 



detS = 1- 



+ 



c 

e2xp 
A^ 



(trn(t) +trnt(i)) + e^^'' (^1 + ^|trn(t)p^ 



\2 , , 

(detn(t) +dettf(i)) - -zg-e^^^ (trnt(t)detn(i) +trn(t)detn+(t)J 



+_e4xp|detn(i)|2, 
2iT= \ 2(|trn(t)|2 - l)e2xp + 2e2x^ 



(44) 



1 + detn(t) + - ^ ) detnt(t) 



X2 

-2— e^^^ 
C 



1 + trnt(t)detn(t) + (^1 - ^ ) trn(t)detn^(t) 
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+ -2(C + A)e^^ + 2^e2>^^trnt(t) + 2^ (^1-^^ e^^^^det nt(t) U{t) 

-2(C - A) - 2^ e^^^tr U{t) + 2^(^1 + e^^^det U{t) 

where the coordinate of the envelope soUton is given by 

Xpix,t) = -2XoC{x-2kt), (46) 

and the time dependence of the phase modulation is embedded in the polarization matrix, 
namely, 

n(t) = ne^^^oCAt_ (47) 

When we take the limit Aq — >■ with AqA and AqC kept finite in eqs. (43)-(45), Q converges 
to the form of eqs. (32)-(35), accompanying the parameters kR = — 2AoA, kj = k and e = 
— ln(4AoA). Here, kn and kj are independent free parameters, apart from the trivial initial 
displacement e. In this sense, in spite of the reduction, we can still regard the form of one- 
soliton solutions (43)- (45) as a general form of those under VBC. 

We can also take another limit. That is, we consider the reduction to the single-component 
case. If we set 

*.o. n^(^^ J. (48) 

the (l,l)-component of Q becomes 

n p-i(2Ai*+5) _ X ox Ccos(4AgCAt + e) + iAsin(4A^CAt + 0) 

Qn-e >-Xo- 2AoC ^ ,„3h(2AoCx + ^) - cos(4AgCAt + 9) ' ^^^^ 

where e^ = C/A. The form (49) was given in ref. 18. We thus verify that our soliton solutions 
are also generalization of those for the single-component NLS equation under NVBC.-*^^^ 
4-3 Spin states 

In this subsection, we discuss the spin states of one-soliton solutions with a finite back- 
ground, by calculating the total spin. The conservation law guarantees that we obtain the 
total spin Ft from integrating at arbitrary time. Therefore, we can select the time so that the 
calculation becomes easier. We concentrate on the case of purely imaginary discrete eigenval- 
ues. As a result, we see that the DW-type is associated with the ferromagnetic state, whereas 
the PS-type is associated with the polar state. One-soliton solutions for purely imaginary dis- 
crete eigenvalues (43)-(45) include those under VBC apart from an initial displacement, and 
therefore the classification about the spin states presented below is wider than that performed 
before.^°'^^) 

4-3.1 Ferromagnetic state 

For detn = (DW-type), we substitute eqs. (43)-(45) into eq. (21) and calculate the total 
spin. The time t' such that trn(t') -l-trn^(t') = is suitable for the calculation. The result is 
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Ft = 4Aor ^ 



/ 2ARe{a*(/3 + 7)} \ 



(50) 



(51) 



-2CIm{a*(/3-7)} 

V - ItP) J 

with the modulus 

|FtP = (4AAo)V. 

The total number of the particles transformed into a soliton, iVr, is calculated by eq. (20), 

iVT = 4AoC, (52) 

and the range of the value taken by l-FxP is expressed in terms of iVr, 

N^<\FTf <N^ + {AXof. (53) 

Remark that, in the vanishing limit Aq — ^ 0, the modulus of the total spin is always equal to 
the total number of the particles, namely, I-FtI ^T- 

With nonzero total spin, the DW-type of solitons belongs to the ferromagnetic state. 
Since inter-atomic ferromagnetic interactions are supposed here, solitons tend to take the 
ferromagnetic state or DW-type. In various contexts of physics, the domain-walls are topo- 
logical solitons related with the symmetry breaking. Here, resulting from the domain-walls, 
the magnetic entity emerges as the spontaneously broken symmetry. It is worthy to notice 
that the case jiYxl < |Ft| may happen. The background is spinlcss, but its internal spin state 
appears to be affected on the ground that the ferromagnetic soliton runs over the background. 
Thereby, the background contributes to the total spin. 

4.3.2 Polar state 

If detn / (PS-type), the solitons show the other magnetic property. The time t' such 
that det n(t') = det 11^ (t') > is suitable for the analysis. After lengthy calculations, the local 
spin density at such time is derived as 

( ) = 8Xle-'^E-\xP') ( ^'"^ ) {ce^"S(xpOsinh(xpO 

-e^ [(C^ - A2)trn(t') + (C^ + A2)trnt(t')] sinh(2xp0 

+ [aVC • ((tr nt(t'))' - |tr n(i')P) + 4Cdetn(iO -I- 2C(|tr n(t')|' - 1)] sinh(xp^ 
+h.c., (54) 
fy = 32AgAe-=^nmK(/? - ^)}E-\xp') 



2C e*' sinh(2xp') - (trn(t') + tvU\t')) sinh(xp' 



(55) 
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where H(xp') is an even function of xp'-, XP'i^^t) is a parallel-shifted coordinate, and is a 
constant: 

EixP') = 2cosh(2xP') - ^(trn(t') +trnt(t'))cosh(xP') 

e + ltvjW 2 
+ A2detn(t') ^^^> 

XP' = XP + v = -2AoC(a: - 2kt') + v, (57) 
f = In 



^(detn(O)^/^ 



(58) 

Note that fx and fz share the same functional form. Each component of the above local spin 
density is an odd function of xp' and, in particular, it has the node at the same point xq such 
that XP'i^Ojt') = 0, namely, xq = 2kt' + (2AoC)~^'*-'- Consequently, the total spin amounts to 
zero: 

Ft = J dxf{xP') = {0,0,Of. (59) 
For this reason, the PS-type of solitons, on the average, belongs to the polar state. 

5. Two-Soliton Collision 

We proceed to the discussion of two-soliton collisions in the integrable spinor model (9). 
Two-soliton solutions are obtained by setting N = 4 (M = 2) in the formula (13). There 
exist two independent discrete eigenvalues and symmetric polarization matrices, respectively, 
i.e., Ai = A2 and Hi = for one of solitons, A3 = A4 and = Il\ for the other. Each 
soliton is separated at t = ±00. Then, a two-soliton is asymptotically two one-solitons. The 
classification of one-soliton solutions based on the values of the determinants of polarization 
matrices, discussed in the previous section, is thus valid for two-soliton solutions. 

The derivation of the explicit form is more complicated than that in the case for one-soliton 
solutions. For the derivative NLS equation under NVBC, explicit two-soliton solutions and 
shifts of soliton positions due to collisions between solitons have been analytically obtained, 
in the case of purely imaginary eigenvalues, where complexity of calculation is considerably 
reduced.^°) This strategy, however, does not stand in our NLS equation under NVBC. The 
reason is understood from eq. (46). In the spinor model, purely imaginary eigenvalues give 
two solitons with the same velocity 2k, and they do not collide with each other. Accordingly, 
we can not investigate the properties of collisions for purely imaginary discrete eigenvalues. 
No one has studied explicit multi-soliton solutions of the NLS equation under NVBC, even in 
the single-component case, because of the computational complexity. 

Here, we graphically show the characteristic behaviors of two-soliton collisions in the 
spinor model, by use of the exact solutions given by the ISM. Referring to them, we carry 
out the qualitative discussions. Although the presented graphs are depicted for the specific 
parameters, much the same behaviors are observed for arbitrarily selected parameters. 
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Figure 2 illustrates the behavior of a mutual eoUision between two PS-types, where 
detlli 7^ and detlla ^ 0. One can see that, in all three components. Both solitons re- 
tain their shapes before and after the collision, which is the common property with solitons 
in the single-component case. In this sense, PS-PS soliton collision is essentially equivalent to 
two-soliton collision of the single-component NLS equation. 

Figure 3 illustrates the behavior of a mutual collision between DW-type and PS-type, 
where det Hi = and det lis ^ 0. The behavior of collisions between DW-type and PS-type is 
qualitatively alien from that between two PS-types. One observes that, in PS-type, much of 
the population initially inhabiting the hyperfine substate \F = l^mp = ±1) is transferred into 
the hyperfine substate \F = \,mp = 0) due to the collision. In contrast, in DW-type, such 
spin transfer does not occur, and the domain-wall shape is preserved against the collision. This 
phenomenon can be interpreted as follows. DW-type, with nonzero spin, can affect the internal 
spin state of PS-type, whereas PS-type, which is expected to have zero spin in total, does not 
affect the internal spin state of DW-type. This kind of spin-transfer phenomenon, called the 
spin- switching, has been first predicted for the case of VBC."*^^^ Due to the conservation laws, 
the total number, the total spin and so on are invariant throughout the collision process. 
Population-mixing among internal degrees of freedom is permitted, as far as the conservation 
laws are not violated. The spin-switching is one of the dynamical processes which make the 
spinor solitons more interesting. 

Finally, for a mutual collision between two DW-types, where det Hi = det lis = 0, the 
shapes of both solitons arc expected to be deformed by the collision since each soliton carries 
nonzero total spin. In fact, drastic population-mixing is seen in Fig. 4, which shows an example 
of this kind of collisions. One finds that domain- walls "repel" at the collision region, rather 
than collide. 



Fig. 2. Density plots of j^i p (a), \(f)o'\^ (b) and P (c) for a mutual collision between two PS-types. 




(a) 



(b) 




The parameters used here are fc = 1, Aq = 1, Ai = 1.03i, A3 = 1.05 -|- i. Hi 





. The velocity of the right (left) moving soliton is 2.00 (—3.41). The collision 
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(b) (c) 
Fig. 3. Density plots of \4'i\'^ (a), |(^oP (b) and (c) for a mutual collision between DW- 

type and PS-type. The parameters used here are the same as those of Fig. 2, except for 
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Fig. 4. Density plots of (a), \(f)o\'^ (b) and (c) for a mutual collision between two DW- 

'1/2 i/2 \ 
i/2 -1/2/ 

1 o' 



types. The parameters used here are the same as those of Fig. 2, except for IIi 







. The values more than 2 are colored white. 



6. Concluding Remarks 

In this paper, we have investigated dynamical properties of matter-wave bright solitons 
with a finite background in the F = 1 spinor Bose-Einstein condensate. To perform our anal- 
ysis concretely, we have exploited an integrable spinor model with a self-focusing nonlinearity 
and the inverse scattering method under nonvanishing boundary conditions. The situation 
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that matter- wave solitons are located on a finite background fits to the experiments. 

One-soliton solutions are derived explicitly and studied in detail. From the point of the 
mathematical view, they offer general forms of bright soliton solutions of the NLS equation. 
We have confirmed that our one-soliton solutions include those obtained in the previous works. 
One-soliton solutions are classified into two kinds by the difference of boundary conditions; 
DW-type and PS-type. The spin density profiles of one-solitons vary depending on the bound- 
ary conditions. In the case of purely imaginary discrete eigenvalues, we have analytically shown 
that DW-type is in the ferromagnetic state, whereas PS-type is in the polar state. The exis- 
tence of two distinct magnetic properties for one-soliton solutions also gives rise to fascinating 
phenomena in the case for two-soliton collisions, for example, the spin-switching. The above 
results for bright solitons with a finite background agree with those for bright solitons under 
VBC^O'ii) and dark solitons.^^) 

Several problems still remain. It is desirable to extend our analysis to the case of general 
discrete eigenvalues. The computations of the conserved quantities other than the total spin 
are also required. (One approach is given in Appendix.) In addition, we wish to investigate 
analytical properties of general A^-soliton solutions under NVBC in the spinor model. Needless 
to say, too complicated calculations are inevitable for the above problems. The remaining 
problems should be discussed elsewhere. 

We conclude that the properties of the multiple matter-wave solitons in the spinor BECs 
are interesting and should be useful in various applications. Bright solitons arc preferable to 
dark solitons for applications, because of the advantage in the propagation distance. We hope 
that our analysis contributes to illuminating dynamical properties of solitons in the spinor 
BECs, which should be demonstrated experimentally in near future. 
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Appendix: Several Conserved Quantities of One-Soliton Solutions 

The conserved quantities help us to understand the dynamics of the system. In this ap- 
pendix, we calculate the total number, the total spin, the total momentum and the total 
energy of the spinor model. We assume that, in addition to purely imaginary discrete eigen- 
values, n is a real symmetric 2x2 matrix. The condition that 11 is a real symmetric matrix 
is inherent in the self-defocusing case, i.e., dark solitons. 

For n = nt, one-soliton solutions of purely imaginary discrete eigenvalues (43)- (45) be- 
come the following form at t = t' = (4n — l)7r/8Ag^A for n = 0, ±1, . . . : 




(A-l) 
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where e'''/^ = X/C,. This form is suitable for calculations, since the imaginary part is separated 
from the real one. One can see clearly that the one-soliton solutions under VBC (32) are located 
on a finite background in the form (A-1). Note that the domain- wall shape is lost even for 
det n = there. 

Several conserved quantities of the solitons (A-1) are calculated by use of eqs. (20)-(23). 
The results for det 11 = are given by 



ATt = 4AoC, 



/ 2a(/3 + 7) \ 





'2 „,2 



2X1) 



Ft = At 

Et = A'tc|(A;2 
and those for det 11 7^ are given by 

A^T = 8A0C, 
Ft = (0,0,0)^, 
Pt = NThk, 
Et = Ntc I (fc2 



|Ft| = At, 



2 a2 



2A^) - ixie 



(A-2) 

(A-3) 

(A-4) 
(A-5) 

(A-6) 
(A-7) 
(A-8) 

(A-9) 



It is intriguing that, for fixed amplitude and discrete eigenvalue. At, -Pt and Et of the 
PS-type (det n 7^ 0) have just twice values as those of the DW-type (det 11 = 0), respectively. 
This enables us to interpret that the PS-type of solitons is a bound state of the two DW-types 
of solitons. Additionally, for fixed amplitude and total number, the total energy Et of the 
DW-type is lower than that of the PS-type: E^"^ - E^^ = < 0, which suggests that 

the DW-type is physically preferable. This result is consistent with inter-atomic ferromagnetic 
interaction, i.e., C2 < 0. 
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